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A search for the — >■ and — ^ A'~''A'~ decays is performed in ~ 300 pb~^ of pp collisions at y/s = 

7TeV collected by the LHCb experiment at the Large Hadron Collider at CERN. The measured limit on the 
branching fraction of the — decay is B{Bg — > < 1.6 X 10~* at 95% confidence level, while 

that of the B° -5- decay is B{B° -> < 5.1 X IQ-^ at 95% confidence level. 

A combination with the 2010 dataset of LHCb yields a hmit of B{B^^ — 5> fi^fi^) < 1.5 X 10~* at 95 % confidence 
level. 



1. Introduction 

Measurements at low energies may provide interesting indirect constraints on the masses of particles that 
are too heavy to be produced directly. This is particularly true for Flavor Changing Neutral Currents (FCNC) 
processes which are highly suppressed in the Standard Model (SM) and can only occur through higher order 
diagrams. 

The SM prediction for the branching fractions (B) of the FCNC decays — > f^^iJ.^ and — > fi^iJ.^ ^ have 
been computed [T] to be B{B° /x+At") = (3.2 ± 0.2) x IQ-^ and B{B° n+fi-) = (0.10 ± 0.01) x 10"^ 
However New Physics (NP) contributions can significantly enhance these values. For instance, within Minimal 
Supersymmetric extensions of the SM (MSSM), in the large tan/3 approximation '5], B{B'^ — )■ ^'^ n~) is found 
to be proportional to ~ tan^ /?, where tan/3 is the ratio of the vacuum expectation values of the two neutral 
CP-even Higgs fields. Therefore it could be strongly enhanced for large values of tan /3. 

The best published limits from the Tevatron at 95%C.L. are obtained using 6.1 fb-i by the DO collabo- 
ration [3 , 6(5^ < 5.1 X 10-8 and using 2 fb^i by the CDF collaboration 14J, S(B° M'^M") 

< 5.8 X 10"^ and B{B^ 1^'^ 1^^) < 1-8 x 10^^. The CDF collaboration has also presented a prehminary 
result [5 with 6.9 fb^^ that lowers the hmit to B{B^ 1^^ l^-") < 6.0 x 10^^ at 95% C.L. . In the same dataset 
an excess of — > /i^/i^ candidates is also observed compatible with a B{B'^ /i+/i^) — (1.8^q g) x 10^^ and 
with an upper hmit B{B° IJ-'^l-i^) < 4.0 x 10"^ at 95% C.L. The CMS collaboration has recently presented 
a preliminary result ,6J with 1.14 fb^^ and set the hmits B{B'^ IJ'~^IJ-~) < 1-9 x 10"* and B{B^ M^M ) 

< 4.6 X 10-9 at 95% C.L. 

The LHCb collaboration has obtained limits [7] comparable with the Tevatron ones, B{B'^ — > /i+/i^) < 
5.4 X 10-8 and B{B" IJ-^I^^) < 1-5 x 10'^ at 95% C.L. based on 37 pb-^ of luminosity collected in the 2010 
run. This document presents an analysis of the LHCb data taken with an integrated luminosity of 300 pb-^ in 
2011. A combination of this data with the data used for the analysis published on the 2010 data based on 37 
pb-^ is also performed. 



2. The LHCb detector 

The LHCb detector [5] is a single-arm forward spectrometer with an acceptance for charged tracks with 
pseudorapidity 77 of 2 < < 5. It consists of a vertex locator, a warm dipole magnet, a tracking system, two 
RICH detectors, a calorimeter system and a muon system. 

Track momenta are measured to a precision of Sp/p ranging 0.35-0.5 %. The RICH system provides charged 
hadron identification in a momentum range 2-100 GeV/c. Typically, kaon identification efficiencies of over 90% 
can be attained for air ^ K fake rate below 10%. The calorimeter system identifies high transverse energy (Et) 
hadron, electron and photon candidates and provides information for the trigger. The muon system provides 
information for the trigger and muon identification with an efficiency of ^ 95 % for a misidentification rate of 
about 1-2 % for momenta above 10 GeV/c. 



In these proceedings the inclusion of charge-conjugate states is implicit. 
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3. Analysis Strategy 

The analysis presented in this document uses about 300 pb^^ of data collected by LHCb between March and 
June 2011 at ^/s = 7TeV. Assuming the SM branching fractions and a^j^ = 75 ± 14 /ib as measured within 
the LHCb acceptance [5], approximately 3.4 (0.32) M"*"/^" {B^ ^ A'^M ) events are expected to be 

reconstructed and selected in this sample. The general structure of the analysis is unchanged with respect to 
Ref. [3, but has been improved in two main points: the optimization of the multi-variate analysis discriminant 
using geometrical and kinematic variables and the use of the RICH particle identification (PID) in the analysis. 

The selection procedure treats signal and control/normalization channels in the same way in order to minimize 
systematic uncertainties. After the selection, each event is given a probability to be signal or background in 
a two-dimensional space defined by the invariant mass and a Boosted Decision Tree (BDT) from the TMVA 
package [10]. The combination of variables entering the BDT is optimized using Monte Carlo (MC) simulation. 
However the probability for a signal or background event to have a given value of the BDT is obtained from data 
using trigger unbiased — > h^h ^ events for signal and sideband B^^^ — > candidates for background. 

The parameters describing the invariant mass line shape of the signal are extracted from data using control 
samples. The average mass values are obtained from B^ — )■ JsT+tt" and B^ — >■ K^K~ exclusive samples. The 
B^ and B^ mass resolutions are estimated by interpolating the ones obtained with the dimuon resonances 
{J /'(p,ip{2S) and T(15'), T(25') and T(3S')) and cross-checked via a fit to the invariant mass distribution of the 
B^ —>■ h~^h ~ inclusive decays and of the B^ — )■ K^Tr~ exclusive decay. 

The number of expected signal events is obtained by normalizing to channels of known branching fractions, 
B+^ J/ipK+ , B°^ J/ip(j) and B° ^ K+n^ that are selected in a way as similar as possible to the signals in 
order to minimize the systematic uncertainty due to the different phase space. The probability for a background 
event to have a given BDT and invariant mass value is obtained by a fit of the mass distribution of events in 
the mass sidebands, in bins of BDT. The two-dimensional space formed by the invariant mass and BDT is 
binned, and for each bin we compute how many events are observed in data, how many signal events are 
expected for a given BF hypothesis and luminosity, and how many background events are expected for a given 
luminosity. The compatibility of the observed distribution of events in all bins with the expected one for a given 
BF hypothesis is computed using the modified frequentist method method [TT] , which allows to exclude a given 
hypothesis at a given confidence level. In order to avoid unconscious biases, the data in the mass region defined 
by AIbo — 60 MeV/c^ and M^o -I- 60 MeV/c^ have been blinded until the completion of the analysis. 

4. Event selection 

The selection aims at reducing the data size to a manageable level while keeping the efficiency on the signal 
as high as possible. The selection consists of loose requirements on track separation from the interaction point, 
decay vertex quality, compatibility of the reconstructed origin of the B meson with the interaction point and 
muon identification. Events passing the selection are considered B'^^^ /^^M^ candidates if their invariant mass 

lies within 60 MeV/c^ of the nominal iJ^*^^ mass. 

A similar selection is applied to the normalization channels, in order to minimize systematic errors in the 
ratio of efficiencies. The B^ h~^h ~ inclusive sample is selected exactly in the same way as the — >■ /J."*"/^" 
signals (apart from the muon identification requirement) in order to be used as the main control sample in the 
calibration of the BDT and the invariant mass. The J/ip -> /i/i decay in the — > J/ipK^ and 5° — ^ J/ '^4' 
normalization channels is selected exactly in the same way as S^^-j — > I-l^ ii'' signals. candidates are required 
to pass some track quality and impact parameter cuts. 

Fiducial cuts are applied to reject fake signal candidates. All tracks from selected candidates are required 
to have a momentum less than ITeV/c. Furthermore, only B candidates with lifetimes less than 5 x tb^^-^ 
are accepted for further analysis. Dimuon candidates coming from elastic di-photon production are effectively 
removed by requiring a minimal transverse momentum of the B candidate of 500 MeV/c. 

5. Signal and background likelihoods 

The discrimination of the signal from the background is achieved through the combination of two independent 
variables: the BDT and the invariant mass. The invariant mass in the search regions (±60 MeV/c^ around the 
B^g-^ masses) is divided into six equal-width bins, and the BDT into four equal-width bins distributed between 
zero and one. 
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Figure 1: Fits to the dimuon invariant mass spectrum in the [2.9, 3.9] GeV/c? range (left) and in [9, 11] GeV/c^ (right). 

The BDT combines information related with the topology and kinematics of the event: the B^^-^ lifetime, 
the minimum impact parameter of the two muons, the distance of closest approach of the two tracks, the i?^^^ 
impact parameter and pt and the isolation of the muons with respect to the other tracks of the event, the 
cosine of the polarization angle, the B isolation, the minimum of the two muons. Those variables have been 
chosen on the basis of the performance on MC samples. The BDT output is found to be fully independent 
of the invariant mass both for signal and for background. The BDT is trained using Monte Carlo simulations 
[B^^^ pi-^H^ for signals and bb — iifiX for background) and calibrated with data using B'^^^ h^h ^ selected 
as the signal events and triggered independently on the signal in order to avoid trigger bias. 



5.1. Calibration of the invariant mass liltelihood for signal 

The invariant mass shape for the signal is parametrized as a Crystal Ball function. The mass average value 
and the mass resolution are calibrated from data and the transition point of the radiative tail is taken from 
Monte Carlo simulation. 

The average mass values are extracted from K^ir^ and i?" K^K^ exclusive candidates and the 

results are the following: 

m(S°) = (5358.0 ± 1.0) MeV/c^ 
m(B°) = (5272.0 ± 1.0) MeV/c^. 

The mass resolutions are extracted from data by a linear interpolation between the measured resolution of 
charmonium and bottomonium resonances decaying to two muons {J/^, "0(28) and T(IS'), T(2S'), T(3S')) 
as shown in Fig. [T] This result has been cross-checked with the mass resolution obtained via full fit to the 
B^ — h^h ~ inclusive decay and the B'^ — >■ K^'k~ exclusive decay. The interpolation yields (jiB^s) — (24.6 ± 

0.2,tat±1.0syst)MeV/c2. 

5.2. Calibration of the BDT likelihood for signal 

The BDT shape for signal events is calibrated using a — > ~ inclusive sample composed only of events 
Triggered Independently of the Signal (TIS)^ to avoid biases introduced by the hadronic trigger. The number of 
-B^ h'^h~ signal events in each BDT bin is extracted by fitting the hh invariant mass distribution in the tttt 
mass hypothesis [121 . Figure [2] depicts fits to the mass distribution of the full sample and for the three highest 
BDT bins for — !• h'^h ^ TIS events is shown. Those yields are used to compute fractional yield of events 



^Here we consider events that are TIS with respect to the two first trigger levels (LO and HLTl) but accept all HLT2 triggers. 
This is due to the lack of statistics for full (LO X HLTl X HLT2) TIS events. Corrections are computed from Monte Carlo 
simulations and cross-checked on control channels. 
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Figure 3: Probability of signal events in bins of BDT obtained from the inclusive sample of TIS h^h ~ events 

(solid squares). The background probability (open circles) is obtained from the events in the sidebands. 



as a function of the BDT bin as shown on Fig. [3] Cross-checks consist of one-dimensional fits to the inclusive 
sample and of fits to B'^^^ — ^ KK, Kit, tttt exclusive samples selected using the K — it separation capability of 
the RICH system. The maximum spread among the different methods has been used as a systematic error in 
the fractional yield determination. 



5.3. Calibration of the combinatorial background Likelihood 

The BDT and invariant mass shapes for the combinatorial background inside the signal regions are calibrated 
with data by interpolating the number of expected events using fit to the invariant mass sidebands in each bin 
of BDT. The lower sideband is potentially polluted by two contributions other than combinatorial background 
that can affect the fit: cascading 6 — > c/i ^iiX decays below 4900MeV/c^ and _B" h^h ^ events with the 
two hadrons misidentified as muons above 5000MeV/c^. Thus the model used to extract the expected number 
of background events is chosen to be a single exponential fitted to the events in the region [4900, 5000] MeV/c^ 
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Figure 4: Distribution of the /i/x invariant mass for events in each BDT bin. The curves show the models used to fit 
the sidebands and extract the expected number of background events in the signal regions. Red: linear fit in [4672, 52f 2] 
MeV/c^ and [54f 8,6000] MeV/c^ sidebands. Blue: single exponential fit in the same ranges as the red curve. Green: single 
exponential fit in the region [4900,5000] MeV/c^ and [54i8,6000] MeV/c^. Magenta: single exponential fit using only the 
upper sideband ([54i8,6000] MeV/c^). The green curve defines the model and ranges used in the analysis. The and 
5° search windows, delimited by the vertical green and orange lines, have not been considered in the fitting procedure. 



and in the upper sideband [5418, 6000] MeV/c^. 

In Fig. [4] the invariant mass distribution for events in each BDT bin is shown and different fit models used 
for estimating the expected number of background events in the signal regions are superimposed. The fraction 
of combinatorial background expected per BDT bins is depicts on Fig. [3j 

5.4. Peaking background Likelihood 

The peaking background coming from h'^h ^ events for vifhich the two hadrons in the final states are 

misidentified as muons has been evaluated from data to be A'^^o = (0.5 ±0.4) events in the mass window and 
Nbo ~ (2.5 ±2.0) events in the B^ one in the whole BDT range. The mass lineshape of the peaking background 
is obtained from a MC sample of double misidentified B^ h'^h~ events and normalized to the number of 
events expected in the two search windows from data. 



6. Normalization 

The number of expected signal events is evaluated by normalizing with three complementary channels of 
known BF, B^ J/ipK^, B^ — J/ipfj} and i?" — > K^tt~. The first two decays have similar trigger and muon 
identification efficiency to the signal but a different number of final-state particles, while the third channel has 
the same two-body topology but cannot be efficiently selected with the muon triggers. 

The number of observed events A'^B0_5.^+p- is translated into a BF via the relation: 

^REC SELIREC TRIGISEL . ^ _ 

f7 '^norm'^norm tnorm J norm j.^ 

O - "norm X SEL|REC TRIG|SEL ^ 'JZT ^ N ~ ^ ^^-8"^^+^- > 

sig sig sig 1 
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where /so and /norm are the probabihties that a 6-quark fragments into a and into the 5-hadron relevant for 
the chosen cahbration mode. We use LHCb measurement for fd/fs = 3.745 ±0.295 [13]. Bnorm is the branching 
fraction and Naom the number of selected events of the normalization channel. The efficiency is separated into 
three factors: e^^^ is the efficiency to reconstruct all the tracks of the decay including the geometrical acceptance 
of the detector; gSEL|REC jg ^j^g efficiency to select the events which have been reconstructed; gTRiG|SEL 
the efficiency of the trigger on reconstructed and selected events. The ratios of reconstruction and selection 
efficiencies are estimated from simulations and checked with data, while the ratios of trigger efficiencies on 
selected events are determined from data |14l . Finally, anorm is the normalization factor. 

Each A^norm IS obtained with a fit to the invariant mass distribution of the normalization channel. The 
invariant mass distributions for 5+ — !• J/ipK^ and 5° — )■ J/iJ^ are shown in Fig. [H] while the — > K'^Tr~ 
yield is obtained via the — > h^h ^ full fit as pictured in Fig. 2 top left. The normalization factors calculated 
using the three channels give compatible results; the final normalization factors are weighted averages and read 
a^o^^^^^ = (9.84 ±0.91) x 10-^° and aso^^^ = (2.59 ±0.15) x 10^1°. 



7. Results 

The expected number of combinatorial background events, the number of peaking background events as well 
as the number of expected signal events assuming a given BF hypothesis are computed for each 24 bins. The 
distribution of the invariant mass in the four BDT bins is shown in Fig. [6] and [t] for — > /i^/i^ candidates and 
B'^ — >■ /i^/^^ candidates, respectively. The observed events are in good agreement with background expectations. 
For B'^ — )■ ^Ji^ the observed events are also in good agreement with the background expectations and the 
presence of B^ — > events according to SM predictions. 

The compatibility of the distribution of events in the invariant mass vs BDT plane with a given branching 
fraction hypothesis is evaluated using the CLg method [TT]. This provides two estimators: CLg is a measure 
of the compatibility of the observed distribution with the signal hypothesis, while CLb is a measure of the 
compatibility with the background-only hypothesis. 

For the B^ M^A*" decay, the distributions of expected CLg values are shown as dashed (black) lines in 
Fig. |8] (left) under the hypothesis to observe background plus SM events. The green shaded areas cover the 
region of ±1(t of compatible observations. The observed CLg as a function of the assumed branching fraction 
is shown as dotted (blue) line. For the B^ — > /i^/i^ decay, the distribution of possible values of CLg is shown 
as dashed (black) line in Fig. [s] (right) under the hypothesis to observe background events only. The observed 
CLg as a function of the assumed branching fraction is shown as a dotted (blue) line. 

The upper limits read: 

i3(B° ^ < 1.3 (1.6) x lO"^ at 90% (95%) c.l., 

B{B° fi+fi-) < 4.2 (5.2) X 10"^ at 90% (95%) C.L., 

while the expected values of the limits are S(i?° -> M^M ) < 1-2 (1-5) x 10~^ under the hypothesis to observe 
background plus SM events and B{B" < 2.4 (3.1) x 10"^ at 90% (95%) C.L. under the hypothesis 

to observe background events only. 
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Figure 6: 5° — > • distribution of selected dimuon events in the invariant mass plane for the four BDT bins. 

The black dots are data, the light blue histogram shows the contribution of the combinatorial background, the green 
histogram shows the contribution of the -B, — >■ h^h~ background and the red filled histogram the contribution of 
-Ba ~^ M^A*" signal events according to the SM rate. The uncertainty on the expectations is shown by the hatched area. 
The uncertainty on data in the first BDT bin is smaller than the size of the dots. 
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Figure 7: — )■ /i^/^ : distribution of selected dimuon events in the invariant mass plane for the four BDT bins. The 
black dots are data, the light blue histogram shows the contribution of the combinatorial background and the green 
histogram shows the contribution of the B'^ — >■ K^h ^ background. The uncertainty on the expectations is shown by the 
hatched area. The uncertainty on data in the first BDT bin is smaller than the size of the dots. 
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Figure 8: Expected distributions of possible values of CLs (dashed black line) under the hypothesis to observe a com- 
bination of background and signal events according to the SM rate in the case of _B° — > ^~ (left) and under the 
hypothesis to observe background only in the case of — >■ fi'^ /J." (right). On both graphs, the green areas cover the 
region of ilcr of compatible observations. The observed distributions of CLs as a function of the assumed branching 
fraction are shown as dotted blue lines for — >■ fi'^ (left) and — >■ fi'^ (right). The measured upper limits at 
90% and 95% C.L. are also shown. 

The comparison of the observed distribution of events vi^ith the expected background distribution resuhs in 
a p-value (l-CU) of 20% (21%) for the 5° ^ ^i+fi- (5° ^ M+M") decays. In the case of B" fi+^i~, the 
slightly low p-value is due to an excess of the observed events in the most sensitive BDT bin with respect to the 
background expectations. A larger data sample will allow to clarify the situation. In the case of — > n~ , 
when a signal is included at the level expected in the Standard Model the p-value increases to 50%. 

A limit on B{B'^ — >■ fX^ fjL~) has been extracted after combining the analysis on 300 pb~^ of 2011 data with 
the published analysis [7] on 37pb~^ of 2010 data; the upper limits improve to: 

6(5° ^ ^+/i-)(2010 + 2011) < 1.2(1.5) X 10"^ at 90 %(95 %) C.L. 

The above 95% C.L. upper limits are still about 4.7 times above the SM predictions for the and 51 times 
for the 



References 

1 A.J. Buras, |arXiv:1 012.144 7; E. Gamiz et al, Phys. Rev. D80 (2009) 014503; A.J. Buras, Phys. Lett. B566 
(2003) 115. 

2 H. Chao-Shang, L. Wei and Y. Qi-Shu, Phys. Rev. D59 (1999) 011701; C. Hamzaoui, M. Pospelov and M. 
Toharia, Phys.Rev.B59 (1999) 095005; S. Rai Choudhury and N. Gaur, Phys. Lett.B451 (1999) 86; K.S. 
Babu and C.F. Kolda, Phys. Rev. Lett.84 (2000) 228. 

3 V. Abazov et al. [DO Collaboration], Phys. Lett.B693 (2010) 539. 

4 T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett.lOO (2008) 101802. 

5 T. Aaltonen et al. [CDF Collaboration], arXiv:1107.2304 (2011). 

6 S. Chatrchyan et al. [CMS Collaboration], arXiv: 1107.5834 (2011). 

7 R. Aaij et al. [LHCb Collaboration], Phys. Lett. B699 (2011), 330. 

8 A. A. Alves et al. [LHCb Collaboration], JINST3 (2008) S08005. 

9 R. Aaij et al. [LHCb Collaboration], Phys. Lett. B 694 (2010) 209. 

10 TMVA, Toolkit for Multi Variate analysis with ROOT, http://tmva.sourceforge.net. 

11 A.L. Read, J. Phys. G28 (2002) 2693; T. Junk, Nucl. Instrum. Methods A434 (1999) 435. 

12 A. Carbone et al.., LHCb-PUB-2009-031. 

13 R. Aaij et al. [LHCb CoUaboration] , LHCb-CONF-201 1-034, arXiv: 1106.4435] 

14 E. Lopez Asamar et al. [LHCb Collaboration], LHCb-PUB-2007-073. 



